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'. Abstract 

The generalized wordlength pattern (GWLP) introduced by Xu and Wu (2001) for an 
arbitrary fractional factorial design allows one to extend the use of the minimum aber- 
ration criterion to such designs. Ai and Zhang (2004) defined the J-characteristics of a 
design and showed that they uniquely determine the design. While both the GWLP and 
the J-characteristics require indexing the levels of each factor by a cyclic group, we see 
t that the definitions carry over with appropriate changes if instead one uses an arbitrary 

abelian group. This means that the original definitions rest on an arbitrary choice of group 
structure. We show that the GWLP of a design is independent of this choice, but that the 
■ J-characteristics are not. We briefly discuss some implications of these results. 

in 
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1 Introduction 

On 



O 



In a regular fractional factorial design D, the quantities 



^ , Ai{D) = the number of defining words of length i 

U ; 

contain useful information about the design. In particular, the smallest index i for which 
Ai(D) > is the resolution of the design. Moreover, one way of comparing two designs having 
k factors and equal resolution is to compare their wordlength patterns (A\, A2, ■ ■ ■ , Af.) [6j [7]. 
The better design is said to have less aberration. 

While nonregular designs no longer have defining words as such, a generalized wordlength 
pattern (GWLP) can be defined for them combinatorially. This was done for two- level designs 
by Tang and Deng [13], and was generalized to arbitrary (possibly mixed-level) designs by Xu 
and Wu [15] using group characters. 

An intermediate computation in the two-level case gives a set of values that Tang and Deng 
called J-characteristics (first introduced in [4J), and Tang [T2] showed that these numbers 
completely determine the design D, somewhat analogous to the way that a defining subgroup 
determines a regular design. Ai and Zhang [1] generalized this to arbitrary designs by looking 
closely at the corresponding computation in [15] , 
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In defining generalized wordlength patterns of arbitrary designs, Xu and Wu [15] assigned 
to the ith factor the cyclic group Z Si , where s, = the number of levels of the factor. While this 
choice is a computational convenience, it is also arbitrary, and in fact the calculation of the 
GWLP can be carried through using other abelian groups as well, as we indicate below. 

This, however, raises the following question for non-prime Sj. Since the (irreducible) char- 
acters of two groups of equal order will generally be different, does the choice of group affect 
either the J-characteristics or the GWLP of a given design? Certainly any dependence of the 
GWLP on an arbitrary choice would raise a serious question about its use in comparing designs 
using relative aberration. It will be clearly seen that the J-characteristics do depend on this 
choice. However, perhaps surprisingly, this does not affect the values of the GWLP. That is 
our main result. 

There are many excellent expositions of character theory, such as [8], [9] and [10] . In 
general we will mention known results without citation. We will also use a number of results 
from multilinear algebra (the theory of tensor products). These are collected in an appendix. 

Notation. We will denote the integers by Z, and the integers modulo s by Z s as above. 
The complex numbers will be denoted by C and complex Euclidean space by C s . Vectors in 
C s will be viewed as columns. The conjugate of z £ C will be denoted by z, the transpose of a 
vector or matrix by a prime ('), and the adjoint (or conjugate transpose) of a matrix or linear 
transformation A by A*. The inner product of v = [v±, . . . , v s ]' and w = [u>i, . . . , w s }' £ C s is 
given by 

s 

(V,W) = ^ViWi. (1) 
i=l 

The cardinality of a set E will be written \E\. 

The Hamming weight of u = (u\, . . . ,Uk), wt(u), is the number of nonzero components of 
u. (In Section [2] we will replace "nonzero" by "nonidentity" in order to deal with groups whose 
identity element is not 0.) 

We alert the reader to the fact that we will use G (or Gj) as an index set, with elements g 
or h. Sometimes such sets will be groups, but often they will be viewed just as sets. We will 
try to make absolutely clear from context when a result requires a group structure and when 
it doesn't. 

2 Definitions 

A fractional factorial design on k factors is a multisubset D of a finite Cartesian product 
G = G% X • • • X Gfc, that is, the set G with the element g repeated 0(g) times, 0(g) > 0. 
The set Gi indexes the Sj levels of factor i, and we let s = s\ ■ ■ ■ Sk- We will refer to O as the 
counting or multiplicity function of D. The elements (fe-tuples) of the design are referred to as 
runs, and the number of runs in the design, counting multiplicities, is 

N = \D\ = Y J 0(g). (2) 

96G 
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The design D may also be viewed as an orthogonal array, particularly if its runs are displayed 
in matrix form, say as columns of a k x TV matrix. 

In [15] Xu and Wu denned the generalized wordlength pattern (Ai(D), . . . , Ak{D)) of D as 
follows. If Gi has Sj elements, we take Gi = Z Si , the additive group of integers modulo Sj. This 
makes G an abelian group. To each g E Z s we associate a function Xg '■ Z s ~~ > C such that 

X 9 W=^, (3) 

where £ is a primitive sth root of unity (say £ = e 27rl//s ). For elements = (gi, ■ ■ ■ , gk) and 
/i = (/ii, . . . , h}~) of G = G\ x • • • x Gk, we let 

and defin^ 

X 9 (D) = Y J 0(h) Xg (h). (5) 

/iGG 

Finally, the "generalized wordlengths" are given by 

A 3 (D)=N- 2 \Xg(D)\ 2 for j = 1, . . . , fc, (6) 

wt(a)=j 

where wt(g) is the Hamming weight of g. 

Ai and Zhang pQ note that when si = ■ ■ ■ = s& = 2 the quantities Xg(D) are the J- 
characteristics of Tang and Deng, and rename them so in the general case, with the notation 

Jg(D). 

We now indicate the way in which other groups may be used in (|5|) and ©. 

The functions Xgi are the irreducible characters of the group Z Si , and so the functions 
X 9 are the irreducible characters of G. Among these is Xe = 1, the trivial character of G, 
corresponding to the identity e of G. Something similar holds for abelian groups, in particular 
the indexing of irreducible characters by group elements. 

Specifically, the irreducible characters of an abelian group G are precisely the homomor- 
phisms of G into the multiplicative group C* = C \ 0. The indexing of these characters is based 
on the following result. 

Theorem 2.1. Let Irr(G) denote the set of irreducible characters of the group G. If G is 
abelian, then Irr(G) forms a group under pointwise multiplication, and if G is also finite, then 
G = Irr(G). In particular, the identity element of G corresponds to the trivial character of G. 

The isomorphism is not canonical - and, in particular, not unique - as it depends on the 
representation of an abelian group as a product of cyclic groups (the Fundamental Theorem 
of Abelian Groups), and for cyclic groups on the choice of root of unity in ([3j). (See, e.g., [U 
Theorem 2.4]). We will assume that we have fixed an isomorphism Gi — > Irr(Gj) for each i, 
and thus an indexing of the irreducible characters of Gi by group elements. We will not need 

1 In [15] Xs(^) i s defined as '}2 heD XaW, and it is to be understood that the hth term is repeated the number 
of times h appears in the design |14) . Equation {5), which is essentially the same as that used in [T|, makes this 
explicit. 
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to know the indexing explicitly. The irreducible characters of the direct product G are still 
given by 

J-characteristics and generalized wordlength counts are still defined by ([5]) and ©, respec- 
tively, where we now define the weight of the element g = (g±, . . . , g^) G G to be the number of 
nonidentity components of g. Our main result is this: 

Theorem 2.2. The quantities Aj(D) in (fSjj are independent of the group structure of G. 

The proof of this theorem is given in Section 01 Before considering this, we take a moment 
to study the effect of the choice of group on the J-characteristics of a design. 



3 J-characteristics. The character table. 

We see that the irreducible characters of a finite abelian group G of order s may be written 
X gi , • • • , Xg a > where gi are the elements of G in some order. The values Xgi (dj) form the character 
table of G, the columns of which are mutually orthogonal and of norm ^/s (with respect to the 
inner product {!])). Another way to say this is that the s x s matrix H formed by this table 
has the property that H*H = HH* = si, where H* is the adjoint of H (H is thus a complex 
Hadamard matrix). 

Let G = Gi x • • • x Gk where each Gi is an abelian group, so that G is as well, and assume 
that the elements of G are ordered in some fashion. (Ai and Zhang pQ use a lexicographic or 
Yates order).) If we consider the set of J-characteristics Xg{D) and the counts 0{g) as s x 1 
vectors x an d O indexed by g E G, then ([5j) may be written 

X = HO. (7) 

Multiplying through by H*, we see that H*x = H*HO = sO, so that 

= (l/s)H* X , 

and in particular that the J-characteristics determine the design. This is Theorem 1 of pQ. 

However, in general H depends on the group structure of G, and so from ([7]) or directly 
from © we see that the values of the J-characteristics depend on the choice of group structure. 
This is illustrated with the following example. 



Example 3.1. Consider the 3- factor design 

0000 aaaa bbbb cccc 

D - Oabc Oabc Oabc Oabc . (8) 

Oabc bOca acOb cbaO 



Each factor has 4 levels, namely 0, a, b, and c, and each column is a treatment combination. 
One can check that this is an orthogonal array of strength 2 and index 1 (it is taken from [5] , 
where it is shown to be non-regular). 

For each factor the symbol set Gi = {0, a, b, c} may be given two group structures, namely 
that of the cyclic group Z4 and that of the "Klein 4- group" V (isomorphic to Z2 x Z2). Table I5TT1 
displays the non-zero values of Xg(D) as g runs over the 64 elements of G = G\ x G 2 x G3, 
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Table 1: J-Characteristics for Design D under two different group structures 
For the design in (JSj) , the value of Xg(D) is given for each g £ G, where G = either 
Z4 x Z4 x Z4 or V x V x V . Those g for which both values of Xg(D) = are omitted. 
(Computation was done in Maple.) 
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cbc 
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6 - 2i 
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-A 


4-4i 
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-4i 
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where the groups G{ are all Z4 or all V. For example, Xaaa(-D) = — 6 — 2i using Z4, but = 8 
using V. Thus we see that the values of the J-characteristics depend on the group structure. 
Note that for both group structures we have Xg(D) = if wt(g) = 1 or 2. Such group elements 
have been omitted from the table for convenience. 



It is not hard to calculate the values Aj(D) given by equation 



(The computation is shortened by the fact that | ±a±6i| 2 
group structures we have A\(D) = A%(D) = and As(D) 



+ b' 



, where we have iV = 16. 
We find that under both 



3, as guaranteed by Theorem 12.21 



Before we leave this topic, we develop the properties of the character table a little further. 

In enumerating the elements of a group G one typically chooses g± = the identity. With 
this convention, which we shall adopt, Xgi 1S the trivial character of G, so that Xgi(h) = 1 for 
all h £ G. On the other hand, since G is abelian, Xg{9i) = 1 f° r every g £ G, and so we see 
that H must have the form 

/ 1 1 ••• 1 \ 

1 * • • • * 

H = 

\ 1 * • • • * J 

The matrix U = (l/y/s)H is said to be the normalized character table of G. We list its 
important properties here, which follow from the preceding. 

Lemma 3.2. U unitary (U*U = UU* = I), and in particular defines an isometry on C s 
((Uv,Uw) = (v,w}). Ife= [1,0,..., 0]' andb= (l/ v /i)[l, ...,!]' then 



Ue 



U*e. 
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4 Independence of group structure 



By imposing a group structure on the set G = G\ x • • • x Gk, we define the irreducible characters 
X g - We want to show that the numbers 

E \Xg(D)\ 2 , j = l,-..,k, 
wt( a )=j 

appearing in ([U]) are independent of the group structure chosen. This sum is somewhat unwieldy, 
and so we will break it into smaller sums over elements g which are not only of weight j but 
also differ from the identity in exactly the same components. 

To begin with, we fix an order of the elements in each set Gi, with the understanding that 
whenever we impose a group structure, the first element will be the identity of the group. We 
may denote by lj the chosen element of Gi. 

Now, for each J C {1, . . . , k} with \ J\ = j, let 

Sj = {g = (gi,--- ,g k )£G: gi ?UiKi<= J}. (9) 

(Here J is merely an index set and has no relation to the J-characteristics mentioned earlier.) 
Clearly the sets Sj are disjoint and their union is the set of elements of G of weight j. Then 

E \x 9 (d)\ 2 =E E\x 9 (D)\ 2 - 

wt(a)=j \J\=j g^Sj 

We will show that for each J the inner sum 

E MD)\ 2 (10) 

is independent of the group structure chosen. 

To do this, we will write these sums as squared norms of elements in an appropriate subspace 
Vj of C s . Assuming a fixed ordering of the elements of G, the components of a vector v G C s 
are complex numbers indexed by the elements of G, something like 

v = [... , Vg ,...]'. (11) 

The standard basis elements are of form 

e g = [0,... ,0,1,0. ..,0]', (12) 

where 1 occurs in just the g-th coordinate. Then 

v = Evv ( 13 ) 

Let 

Vj = { v =[-.., v g ,•••]' G C s : v g = if g i Sj] . 

It is clear that dimVj = \Sj\ = riiej( s « — 1)> an d that the sum ([TDD is II m j(x)H 2 where Mj 
is the orthogonal projection of C s onto Vj. Now the next result follows immediately from (JT]) 
and the fact that H = y/s U. 
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Proposition 4.1. Suppose that G = G\ x • • • x G k where Gi is an abelian group. Then the 
sum in |77J) is equal to 

s\\MjUO\\\ (14) 

where Mj is the orthogonal projection of V s on the sub space Vj, U is the normalized character 
table of G, and O is the vector of multiplicities of the design D. 

Our goal is now to show that the quantity (|14p is independent of the group structure of G. 
A very useful way to describe Vj is as follows. Associate to Gi the Euclidean space C s \ 

(i) 

where the components of a vector v are indexed by the elements of Gi. Let e g be the unit 
vector in C Si having a 1 in the gth place and zeros elsewhere, so that = [1, 0, • • • ,0]'. Define 
the subspaces Vi C C Si , i = 1, . . . , k, by setting 

Vi =(e?) ± , ieJ, 
= span(e^), i $ J, 

where orthocomplement ("*-) and span are within C s \ Thus the vectors of Vi have a zero in the 
first position if i £ J and zeros in all the other positions if i ^ J. 

Let Pi be the projection of C Si onto span(e^), Ii the identity matrix, and Qi = I{ — P%. 

Proposition 4.2. With the above definitions, we have 

Vj = Vi®---®V k . (15) 

The orthogonal projection Mj of C s on Vj is given by 

Mj = Mx®---®M h , (16) 

where Mi is the orthogonal projection of C Si on Vi . We have 

Mi =Qi, i € J, 
= Pi, it J. 

Proof. The vectors in this tensor product are sums of vectors of the form 

v\ <g) • • • (8) Vk, Vi G Vi. 

It is not hard to see that a vector of this form has zeros in exactly the positions indexed by 
g $ Sj, so that V\ <g) • • • (8) V k C Vj. However, 

dim Vi = Si — 1, i £ J, 

= 1 otherwise, 

so dim Vi (8) • • • <g> V k = dimFi • • • dim V k = Y\i. e j( s i ~ !) = dim Vj. Thus ([TSJ holds, and (JTBD 
follows immediately. The formula for Mj is obvious. □ 

We also note the following, which is implicit in equation @. 

Proposition 4.3. If Gi is a finite group having character table Hi and normalized table Ui, 
then G = G\ x ■ ■ ■ x G k has character table H = H\ (8 • • • <8> H k and normalized character table 
U = Ui (8 • • • (8 U k . 
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We use this to evaluate the vector MjUO appearing in (fT4"j) . As in (fT3j) . the vector O of 
multiplicities may be written 

= Y J 0{g)e g . 

geG 

But if g = (g%, . . . ,g k ) € Gi X • • • x G k , then 



where e!- is the unit vector in C Si having a 1 in the jth place and zeros elsewhere. Then 



,0 

id tuc: unit vctuui 



g=(gi,-,9k) 

Thus 

MjC/O = £ 0( 5 ) MjC/(e« ® • • • ® e£ } ) 
g=(gi,-,9k) 

= Y, 0{g)M l U 1 {e$)®---®M k U k {eW). (17) 

9=(9l,—,9k) 

To analyze the (squared) norm of this, we need to analyze the terms in such sums. This leads 
to evaluating MjL/j on the basis elements e 9i . We will just need to do this when Mj = Pi. 

Lemma 4.4. For each i and for every g £ Gi we have 

'Si 



Proof. For simplicity, suppress the index i. Now for any w £ C s we have 

Pw = (w,ei )ei, 

so from Lemma 13.21 we have 

PU v = (Uv,ei)ei = { v, U*e\ )e\ = {v,b)e\ 
for any v, with b = (l/y / s)[l, . . . , 1]'. In particular, 



PU(e g ) = (e g ,b)e 1 = — e 1; 



1 

7s 



as claimed. □ 

We now evaluate the squared norm of sums of form (|17p . This will rest on the following 
calculation. 
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Lemma 4.5. Let Ij be identity matrix of order Sj, and let 

c 9i-9> = E °(^' 

9i+l,---,9k 

the sum of the numbers 0{g) over those g £ G with the first i values fixed at (g\, . . . , gi). Then 
for every < i < k, 

\\(h®---®Ii®P i+1 ®---®P k )UO\\ 2 = V 4.. (18) 

In particular, this quantity is independent of the group structure on G. 
Proof. We see that (I x <g • • • <g Ij <g> P i+ i (8) • • • (g> Pk)UO 

= W e $) ® ■ ■ ■ ® W e £?) ® ® • • • ® PkU k {e^) 

g=(gi,—:9k) 

= J2 °(9) U 1 {e g l ))^---®U i {e g i }) ® — ^e? +1) ® • • • -^ef } 



g={9i,-,9k) 

I £ ( E O( 5 ))C/ 1 (e( 1 i ))®...0C/ J (e«)®ei i+1) ®...®ef ) 

V 1 (91,-,Si) \(ffi+i.—,flfc) / 



{91,-,9i) 



E c si- 9! ^(eff) ® • • • ® ^( e 2) ® e i 1} ® • • • ® 4 \ 



But for each j, the set {Uj(eg), g £ Gj} is orthonormal in C Sj as the unit vectors e^\g G Gj, 
are orthonormal and CTj is an isometry. Hence the elements C/i(e^) 8) • • • (g Ui{e g 1 }) <g e[* +1 ^ 
•••(g) ep are orthonormal in C s , and so ||(Ji <g ■ ■ ■ <g) Jj (g) P i+ i <g • • • (g) P fc )C/0|| 2 



(51,— iff.) 

But this 



1 


Sj+1 ■ ■ 


• s k 


1 




Si+l • • 


■ s k 


1 




• • 


■ s k 



c 2 gi ... gi l|f/i(e«)l| 2 ---|l^(e«)|| 2 [lei^ 1 ^ 

(Sl,-,Si) 



V c 2 

Z_> 9l -9i 

(flli-iffi) 

as all the norms in the next-to-last line are 1. This is formula (|18|) . □ 

Now fix J C {1, ... , k}. 
Proposition 4.6. ||Mji70|| 2 is independent of the group structure of G. 
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Proof. Actually, we will prove something more general, namely that the proposition holds for 
projections M made up of a tensor product of Pj's, Qi's and Ij's, where Qi = h — Pi- Letting 
q = the number of factors Qi in the projection, we prove this by induction on q. 
We simplify matters by proving our result for projections of form 

M = Qi ® • • • Q q ® I q+i (8) • • • (8) J 9 +i (8> Pg+i+i ® • • • <8> P*. (19) 

The proof is the same for projections with other ordering of the tensor factors. 

The base case (q = 0) is precisely Lemma 14,51 For the induction step, assume that the 
result holds for projections having q — 1 factors Q (not necessarily the first q — 1 factors). Now 
Q q = I q — P q , so the projection (fl"9j) is 

M = Ql ® • • • <g> (8) I g <g> I q+1 (8 • • • (8 Ig+i (8 Pq+i+l <8> • • • <8> Pk 

-Ql (8 • • • <8 Qg-l «> P g <8> /g+l (8 • • • (8 Ig+i <8 Pg+i+1 <8 • • • <8 Pfe 

= T 1 -T 2 , 

say. Since Ti = M + T2 and M and T2 are orthogonal, the Pythagorean Theorem gives 

||MC/0|| 2 = \\TtUOf - \\T 2 UO\\ 2 . (20) 

But since T\ and T2 contain q — 1 factors Qj, the induction hypothesis applies to both terms 
on the right-hand-side of (|20f) . and therefore to the left-hand-side, as desired. □ 

By Proposition 14.11 this shows that the sum (1101) . and therefore the quantities Aj(D), are 
independent of the group structure of G. Theorem 12.21 is now proved. 

5 Conclusion 

The definition of the generalized wordlength pattern (GWLP) given in [15] makes sense if one 
chooses abelian rather than cyclic groups to index the levels of each factor. The choice to use 
cyclic groups in [15J is arbitrary, and we have shown that while it does affect the so-called 
J-characteristics of a design, it does not affect the GWLP. This removes a possible ambiguity 
in the definition of the GWLP, and therefore in the use of minimum aberration as an optimality 
criteria for nonregular designs. The choice of cyclic groups may be useful computationally as 
the irreducible characters are then especially simple. 

A special case of the invariance with respect to group structure is already implicit in the 
coding literature [3]. (The connection with regular designs is given in |15j.) However, this 
covers designs in which (a) the index sets Gi are the same (the alphabet) and (b) the design is 
actually a subset of G (so that the counting function O is simply an indicator function). Our 
Theorem 12.21 is quite general, and makes no use of concepts borrowed from coding theory. 

The wordlength pattern of a regular design does not determine the design, and in particular 
does not tell us its alias structure. For that, one needs the defining words. We have seen that an 
analog of the set of defining words of a nonregular design is the set of J-characteristics, at least 
in respect of determining the design. However, as we noted in Section [3j the J-characteristics 
vary with the choice of group structure assigned to factors. Certainly the aliasing structure of a 
design does not depend on this arbitrary choice. The GWLP is independent of this choice, and 
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one may therefore ask just what statistical information it carries. This is a question worthy of 
further investigation. 
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A Multilinear background 

In this section we briefly review some results on tensor products that we have used in this 
paper. We only deal with Euclidean spaces (specifically C k ) since that is all we need here. For 
simplicity we concentrate on the bilinear case (two tensor factors). 

There are many expositions of multilinear algebra, such as that in [2]. An interesting 
exposition with some statistical applications is given in |llj . 

As is well-known, the Kronecker or tensor product of the matrices A (m x n) and B is 



a m iB 



a\ n B 



Om.n. B 



For vectors !)6C° and toeC we thus have 



V (X) w 



V\W 



v a w 



G C ab , 



where v = [v\, ... ,v a ]'. This product satisfies the usual bilinear properties, for example, cv (g) 
w = c{v ® w) = v (g> cw (c a scalar) and A <g) (B + C) = A(&B + A(&C. 

If V C C a and W C C b are subspaces, then we define their tensor product to be the subspace 
of C ab given by 

V <8> W = spanju (g> w : v £ V, w £ W}. 

(Technically, V (8) W is constructed as a free vector space modulo bilinear relations, and is only 
isomorphic to a subspace of C ab , but we will identify it with that subspace.) If {e\, ■ ■ ■ , e^} is 
a basis of V and {/i, • • • , fi} is a basis of W, then 

{ei®fji i = l,...k, j = !,...£} 

is a basis of V <S> W. Thus in particular 

dim(V ® W) = dim(V) • dim(W). 



If we use {v\,v%) to denote the inner (or dot) product and = ( v , v ) the norm, then 
we have 



(Vi <g)Wi, V2®W 2 ) = (vi,V 2 )(wi,W2), 



and in particular 



\v (&) w\ 



\v\\ \\w\ 
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The norm and inner product of vectors of the form ^ vi <8> wi are calculated by expanding the 
inner product in the usual way. 

If Ti is a linear transformation on Vi, then T = T\ % T 2 is a linear transformation on V\ <8> V2 
such that 

T(v 1 ®v 2 )=T 1 (v 1 )®T 2 (v 2 ). 

T is evaluated on sums of such terms by linearity. The matrix of T is given by the Kronecker 
product of the matrices Tj. Finally, if S = Si (g> S2 is a linear transformation such that 5jTj is 
defined for each i, then 

ST = 5iTi 5 2 T 2 . 

All of the preceding extends in the obvious way to more than two tensor factors. 
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